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ABSTRACT 

Fourier-transform, mid-infrared absorption spectra of mixed argon/naphthalene matrices at 5 K are shown with 
ratios of argon-to-naphthalene that vary from 1000 to 0. These spectra show the changes as naphthalene clustering 
and aggregation occurs, with moderate spectral shifts affecting the C-H vibrational modes and relatively small or 
no shifts to the C-C and C-C-C vibrational modes. The possible contribution of homogeneous naphthalene clusters 
to the interstellar unidentified infrared bands is discussed. The contribution of polycyclic aromatic hydrocarbon 
(PAH) clusters to the 7.7 /xm emission plateau and the blue shading of the 12.7 /i m emission band are identified 
as promising candidates for future research. In addition, since PAH clusters are model components of Jupiter and 
Titan’s atmospheres, the information presented here may also be applicable to the spectroscopy of these objects. 
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1. INTRODUCTION 

A common phenomenon associated with the interstellar 
medium is the emission of unidentified infrared bands (or 
UIBs) at 3.3, 6.2, 7.7, 8.6, 11.2, and 16.4 /xm (3030, 1610, 
1300, 1160, 890, and 610 cm -1 , respectively). The emitters of 
these bands are generally associated with gas-phase aromatic 
molecules, to date, primarily neutral and ionized polycyclic 
aromatic hydrocarbons or PAHs (see, for example, Peeters et al. 
2004; Tielens 2008; Ricks et al. 2009; Boersma et al. 2010, 
and references therein) and possibly protonated PAHs (Ricks 
et al. 2009; Knorke et al. 2009) and aromatic heterocycles 
(Hudgins et al. 2005; Mattioda et al. 2003, 2005, 2008). 
Efforts to synthesize the UIBs typically calculate emission 
spectra under interstellar conditions from a large database of 
experimentally measured or theoretically computed infrared 
absorption measurements of aromatic molecules (Peeters et al. 
2004; Mulas et al. 2006; Bauschlicher et al. 2008, 2009). The 
match between synthesized spectra and observations is close 
enough that the UIBs are now commonly referred to as the PAH 
bands or Aromatic Infrared Bands (or AIBs). 

The AIBs overlap a series of wider emission “plateaus” 
which vary independently from the emission bands (Bregman 
et al. 1989; Boersma et al. 2009; Tielens 2008) and which 
are attributed to clusters of PAH molecules (Peeters et al. 
2006; Kimura et al. 2007). The identity of the plateau clusters 
is significantly less well known due to the relative lack of 
experimental data for PAH clusters with respect to individual 
PAHs (see the discussion in Rapacioli et al. 2006 for example), 
the problem of accurately accounting for dispersive interactions 
in quantum chemical models of PAH clusters (Tsuzuki et al. 
2004; Saeki et al. 2006; Rapacioli et al. 2009), and the problem 
of determining accurate nucleation and evaporation rates for 
PAH clusters in interstellar environments (Rapacioli et al. 2006). 
In the case of PAH clusters, this lack of data is particularly acute 
since these clusters are important both as a potential emitter of 
the emission plateaus and as a source of individual PAHs through 
evaporation of their sub-units (Rapacioli et al. 2006; Joblin et al. 
2009). 


Here we report an experimental study on naphthalene clusters 
to start to provide this information. This paper is a first 
step toward addressing the infrared properties of naphthalene 
clusters through transmission measurements of the infrared 
absorption spectra of naphthalene molecules in various degrees 
of aggregation within an argon ice matrix. Naphthalene is, of 
course, a very small PAH by interstellar standards. However, 
this makes it particularly attractive from a theoretical point of 
view since it provides empirical data on a system which is small 
enough to make it computationally tractable. In the absence of 
a complete theory of the nucleation of PAH clusters, it is not yet 
clear how the naphthalene cluster structures in solid solutions 
relate to those in molecular beam experiments and how either 
of these relate to structures of interstellar naphthalene clusters. 
One motivation for this paper is to provide empirical data useful 
for progress in this direction. 

Naphthalene is also distinguished for study by recent ob- 
servational results suggesting that some interstellar environ- 
ments have sufficiently attenuated UV radiation fields that per- 
mit the survival of small PAHs. In particular, Iglesias-Groth 
et al. (2008) claimed detection of the naphthalene cation in the 
Perseus molecular cloud complex while Vijh et al. (2005) at- 
tribute the blue luminescence phenomenon in the Red Rectangle 
to three- and four-ringed neutral PAH molecules. It seems rea- 
sonable that a two-ringed PAH such as naphthalene might also 
be present in a region where three- and four-ringed PAHs are 
indicated. 

Homogeneous naphthalene clusters have been well studied 
using multiphoton ionization spectroscopy (Wessel & Syage 
1990; Benharash et al. 1999; Fujiwara & Lim 2003; Saigusa 
& Lim 1995; Gillieron et al. 2007), mass-selective ionization- 
loss stimulation Raman spectroscopy (Schaeffer et al. 1995; 
Benharash et al. 1999; Kim et al. 1999), and photoelectron 
spectroscopy (Ando et al. 2008) but systematic studies of their 
infrared absorption in the AIB range have yet to be conducted. 
The structures of neutral naphthalene clusters as large as the 
tetramer have been determined for molecular beam experiments 
(see Gillieron et al. 2007 for a recent discussion, and references 
therein for more detail). 
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Table 1 

Summary of Molecular Data for the Spectra Shown in Figures 1-3 


Ar/Nap Ratio 

1000 + 300 

530 ± 90 

200 ± 50 

160 + 40 

47 + 7 

6 + 2 

Argon thickness (fim) 
Naphthalene density 3 
Monomer fraction 13 

6.1 ± 1.5 
2.6 ± 0.3 
0.90 ± 0.03 

10.0 ±0.8 
8.4+ 1.2 
0.82 ± 0.03 

3.3 ± 0.5 
7.2+ 1.4 
0.59 ± 0.08 

7.3 ± 1.7 
20 + 2 
0.51 ±0.09 

2.1 ±0.3 
20+ 1 
0.10 + 0.04 

2.2 ±0.8 
160 + 20 


Notes. 

a Column density times 1 0 H mole cm~ 2 . 
b See the text for the discussion of calculations. 


At present, the best available calculations of the mid-infrared 
vibrational frequencies of naphthalene dimers that we are aware 
of are those of Saeki et al. (2006). Using the GAUSSIAN 03 pro- 
gram at the MP2/cc-pVDZ level, they found four dimer struc- 
tures in the local minimum. Three of the dimer structures are 
those in which the naphthalene molecules are stacked, maximiz- 
ing 7r-7r orbital overlap. These were the most stable. Depending 
on the precise level of theory, these were found to have bind- 
ing energies of 5-10 kcal mol -1 — energies comparable to the 
hydrogen bonds in water ice. A T-shaped dimer structure was 
also computed. This has a considerably weaker binding energy 
which lies between 3.5 and 6 kcal mol -1 . The intramolecular 
frequencies of the naphthalene molecules comprising the four 
dimer structures all show qualitatively similar shifts relative to 
the IR absorption bands. Another motivation for this paper is 
to provide experimental data to refine these quantum chemical 
calculations. 

2. EXPERIMENT 

The infrared spectra presented here were obtained using the 
matrix-isolation technique in which a guest molecular species is 
trapped in a dilute solid solution of an IR-inactive host matrix. 
Here, varying concentrations of naphthalene were trapped 
within an argon matrix condensed upon an IR-transparent 
window for transmittance absorption spectroscopy. 

The substrate is a 0.25 cm thick, 1 inch diameter cesium 
iodide window mounted on the cold head of a Sumitomo model 
RDK-408D cryostat with closed-cycle helium refrigerator. In 
these experiments, the cesium iodide window was maintained at 
5 K throughout. The disk temperature was measured by a gold/ 
chromel thermocouple attached to the oxygen-free high thermal 
conductivity (OFHC) copper mount of the disk. The cold head 
is normally positioned to place the Csl window perpendicular 
to the infrared beam of a Digilab model TFS-4000 Fourier- 
transform infrared (FTIR) spectrometer for infrared absorption 
measurements. The cold head can be rotated 90° to place the 
Csl window facing the argon and naphthalene deposition port. 

The naphthalene sample consisted of roughly 5 g of Aldrich 
99%+ scintillation grade naphthalene kept at room temperature 
and dried by exposure to vacuum. The evaporating naphthalene 
was flowed with Matheson ultra-high purity grade argon sup- 
plied from a separate cylinder through 1 /4 inch stainless steel 
tubing. The mixture was introduced into the vacuum system 
through an Exair model 1010SS cylindrical nozzle aimed at the 
Csl surface at roughly 45° to the surface normal. 

Individual spectra were obtained as follows. An absorption 
spectrum of the 5 K window prior to deposition serves as the 
reference spectrum. The Csl window is then rotated to face 
the gas deposition port for codeposition of naphthalene and 
argon. After the matrix deposition is completed, the Csl window 
is rotated to its spectroscopic position and a post-deposition 


absorption spectmm is measured. The gas flowing from the 
nozzle, once established, was not altered during an experiment in 
order to better produce a uniform deposit. Spectra were obtained 
by co-adding 512 scans with an FTIR resolution of 0.5 cm -1 . 

The argon-to-naphthalene (denoted Ar/Nap hereafter) ratio 
was determined as follows. The argon deposit thickness d was 
estimated from the interference fringes observed in the infrared 
absorption spectra using 

m 

d — — , 
nv 

where v is the frequency of a maximum or minimum in the 
interference fringes and m is an integer for a maximum and a 
half-integer for a minimum. Roux et al. (1980) determined the 
refractive index n of thin solid argon films at 20 K to be 1 .225 
between 700 and 3700 cm -1 . This value was adopted for the 
argon thickness measurements despite a slight change in the 
argon ice density at the 5 K temperature of these experiments. 
The column density of argon is calculated assuming an argon 
density of 1.77 g cm -3 at 5 K (Dobbs and Jones 1957). 

The column density of naphthalene in the argon matrix was 
estimated using the equation 

f r dv 

N = , 

A 

where r is the spectral absorbance of an observed band, v is the 
frequency, and A is the integrated absorbance value calculated 
for the analogous band in fully isolated naphthalene published 
by Hudgins et al. (1994). The average column density derived 
from all bands below ~1500 cm -1 is reported here. Results 
of the argon matrix thickness and naphthalene concentration 
calculations are summarized in Table 1 . 

3. RESULTS AND DISCUSSION 

Figures 1-3 show the vibration-specific regions of the mid- 
infrared spectra of matrices with Ar/Nap ratios of approxi- 
mately 1000, 530, 200, 160, 47, and 6 along with the spectrum of 
pure naphthalene measured using the same apparatus. The spec- 
tra are presented on normalized scales to highlight the spectral 
changes that occur as the naphthalene dilution of the samples is 
increased from highly isolated (Ar/Nap = 1000) to highly ag- 
gregated (Ar/Nap = 6 and pure naphthalene). The positions and 
full widths at half-height (FWHH) of the unsealed absorption 
bands are listed in Table 2 along with peak positions for naph- 
thalene theoretically predicted using Density functional theory 
(DFT). Throughout these discussions, band positions will be 
rounded off to the nearest wavenumber and FWHH to 0. 1 cm -1 . 
Peak positions are given in Table 2 to 0.1 cm -1 . 

Perusal of Figures 1-3 and the band positions in Table 2 show 
that spectral changes are induced by naphthalene clustering. 
The extent of band shift, band broadening, and effect on cross 
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Figure 1. The 3200-2900 cm - 1 absorbance spectra of argon/naphthalene matri- 
ces (Ar/naphthalene ~ 1000, 530, 200, 160, 47, and 6) and the 3200-2900 cm -1 
spectrum of pure naphthalene. Traces are normalized so that the higher fre- 
quency b 2 u band (see the text) has unit height and are vertically shifted for 
clarity. 
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Figure 3. The 1550-900 cm -1 absorbance spectra of argon/naphthalene matri- 
ces (Ar/naphthalene ~1000, 530, 200, 160, 47, and 6) and the 1550-900 cm -1 
spectrum of pure naphthalene. Traces are normalized so that the peak at 
~ 1 392 cm ~ has unit height and are vertically shifted for clarity. 
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Figure 2. The 900-400 cm~ 1 absorbance spectra of argon/naphthalene matrices 
(Ar/naphthalene ~1000. 530, 200, 160, 47, and 6) and the 900-400 cm -1 
spectrum of pure naphthalene. Traces are normalized so that the strongest 
component of the CH in-plane bending complex at ~783 cm -1 has unit height 
and are vertically shifted for clarity. 


section depend somewhat on the molecular vibration involved. 
Both the in-plane and out-of-plane pure C-H vibrational modes 
are more strongly influenced than those involving pure C-C 
or mixed C-C and C-H vibrations. Consequently, the results 
are presented in several subsections based on the nature of the 
affected vibrations. 

The fraction of naphthalene molecules existing as isolated 
molecules within the argon matrix as a function of the naph- 
thalene volume fraction was estimated using a computer model 
based upon the algorithm of Mansson & Rudemo (2002). Each 
naphthalene molecule is modeled as a three-dimensional, oblate, 
hard spherocylinder (or OHSC), which is the set of points that 
are at a distance R from a two-dimensional, planar disk of radius 
a (Cuetos and Martmez-Haya 2008). The van der Waals dimen- 
sions of a naphthalene molecule are given as 8. 12 A x 7.36 A 


x 3.4 A in Hallmark et al. (1991). We choose R — 1.7 A and 
a —3.23 A, corresponding to an OHSC thickness of 2R — 3.4 A 
and outer radius R + a — 4.93 A (the mean caliper radius of an 
8. 12 A x 7.36 A rectangle), to reproduce the van der Waals’ 
thickness and mean radius of a naphthalene molecule. 

The steps of the computation for each simulation are as 
follows. First, a user-input intensity parameter for a Poisson 
distribution was used to generate the random number of trial 
OHSCs within a 200 A cube. The center of each OHSC was 
randomly positioned within this cube and given a random unit 
vector for its orientation in space. Probability distributions for 
the positions and orientations were uniformly distributed over 
the full cube and the full set of unit vectors respectively with 
independent values given to each OHSC. 

Second, the global thinning procedure of Mansson & Rudemo 
(2002) was carried out. Each trial OHSC is given a random, 
independent mark from a uniform probability distribution over 
the range 0-1. If any two OHSCs are determined to be 
intersecting, the one with the higher mark is retained and the 
other is deleted. Both OHSCs are deleted if their marks are equal. 
Distances between OHSCs were computed using the iterative 
procedure described by Cuetos and Martmez-Haya (2008) to a 
tolerance of 10” 4 A. 

To evade edge effects, only the surviving OHSCs within the 
central 140 A cube are used to compute statistics. The volume 
fraction generated by a set of 1000 simulations is given by 
the average number of the surviving OHSCs per unit volume 
inside of this central cube times the OHSC volume. The 
monomer fraction of a set is taken to be the total fraction of 
these survivors of that set that had no nearest neighbor closer 
than 3.76 A (the diameter of an argon atom). The volume fraction 
and the monomer fraction were calculated for 10 sets of 1000 
simulations in order to calculate a mean value and a standard 
error of the mean. These values were determined for a range of 
initial intensity parameters between 1 0 fl and 10 -3 A” 3 . 

The monomer fraction as a function of 1 minus the volume 
fraction for this set of initial intensity parameters is shown as 
the filled squares in the log-log plot in Figure 4. These points 
are well represented by a power law with power 26.3 ± 0.5. 
Since each naphthalene molecule replaces approximately four 
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Table 2 

Experimentally Measured Naphthalene Band Positions (FWHH in Parentheses) Compared to the Theoretically Computed Spectrum Using DFT 


Ar/Nap = 1000 

Ar/Nap = 530 

Ar/Nap = 200 

Ar/Nap =160 

Ar/Nap = 47 

Ar/Nap = 6 

DFT a 



C-C-C and C 

-H out-of-plane bending region 





475.2(1.7) 

479.1 (6.2) 

478.8(8.9) 

478.6 (7.5) 

480.0 


621.1 (4.0) 

619.6(2.2) 

620.1 (3.2) 

619.6(2.2) 

618.9(1.9) 



624.2 (0.5) 

620.8 (1.0) 

620.8 



632.1 

783.6 (0.9) 

783.6 (0.9) 

783.6 (1.3) 

783.6(1.7) 

783.6 (9.2) 


788.4 

787.2 (2.6) 

787.2(5.1) 

787.2(3.0) 

787.2 (7.7) 

789.3 

787.2 

786.2(8.5) 


C-C stretching and C-H in-plane bending region 



960.1 (6.8) 

961.8 (11.2) 

961.3 (3.8) 

963.4 (2.5) 

962.5 (8.0) 

963.5 

1012.6(1.4) 

1011.7 (2.6) 

1012.1 (1.4) 

1012.1 (4.2) b 

1012.1 (2.7) 

1011.9(2.9) 

1009.7 

1016.5 (1.2) 

1016.5 (1.4) 

1016.2(1.4) 

1016.0 (2.2) b 

1016.0 

1051.2(7.4) 

1051.0(12.2) 





1085.2(19.6) 

1078.7(26.2) 

1077.5 (13.4) 


1129.3 (0.7) 

1129.3 (1.0) 

1129.3 (0.1) 

1129.1 (1.4) 

1129.3 (2.7) b 

1130.5 (3.5) 


1131.0(1.0) 

1131.0 (3.0) b 

1131.0 (3.3) b 

1131.3 (2.3) b 

1131.3 (3.2) b 


1131.9 




1140.7(2.2) 

1139.5 (5.9) 

1139.5 (5.0) b 
1148.1 (3.9) 






1161.6(0.6) 

1165.3 (4.1) 

1169.7 



1212.5 (1.2) c 

1212.5 

1212.0 (2.4) b 

1211.8 (3.5) 

1208.8 


1214.4(0.5) 

1214.0(1.7) 

1213.7 (4.1) 

1213.7 (2.6) b 






1245.5 (1.0) 

1243.8 (13.2) 

1243.6(7.8) 


1273.0 (0.9) 

1272.8 (1.4) 

1272.8 (4.2) 

1272.8 (3.3) 

1270.1 (5.0) 
1345.8 (1.3) 

1270.1 (4.4) 

1270.0 





1346.6 

1347.8 (9.1) 



1362.0(1.0) 

1362.2(2.0) 

1362.0(1.7) 

1362.0(2.3) 

1362.7 (3.1) 

1358.3 

1392.1 (1.0) 

1392.1 (1.0) 

1391.9(1.9) 

1391.9(1.7) 

1391.9 (2.6) 

1391.6(3.2) 

1399.8 





1508.6 

1508.1 (4.3) b 

1507.9 

1513.2(0.9) 

1513.2(1.5) 

1512.9(1.7) 

1512.9 (2.4) 

1512.4(2.2) 

1512.4(2.9) 






1525.9(9.4) 

1525.5 (3.6) 




C- 

H stretching region 









2926.7(5.1) 

2986.4 (3.8) 

3009.4 (6.5) 


3029.2 (15.2) b 

3027.8(5.4) 

3029.7 (18.8) d 

3028.7 (3.7) 

3027.8(4.9) 

3026.6 (5.3) 


3042.5 (12.9) c 

3041.0(3.3) 

3040.8 

3041.5 (6.7) 

3041.0(4.0) 

3039.6 (5.0) 

3043.5 

3045.4 



3062.0 





3065.6 (3.8) 

3065.1 (4.1) 

3065.4(7.6) 

3064.9 (10.9) 

3063.0(10.4) 

3059.1 (12.6) 

3062.5 

3078.9 (2.9) 

3079.1 (3.2) 

3078.4 (4.0) 

3077.7 (7.3) 

3075.7(7.1) 

3073.8 (7.2) 
3096.5 (5.7) 
3105.4 (6.9) b 

3076.5 


Notes. Units of cm 1 for all entries. 
a DFT calculations from Langhoff (1996). 

b Peak is only partially resolved. FWHF1 measured up to interpeak minimum. 
c Peak is just barely resolved. FWHFI measured up to interpeak minimum. 
d Quoted FWFIFf is for entire complex of partially resolved peaks. 


argon atoms within the matrix (Crepin et al. 2001), we use 
the naphthalene/argon ratio times four as the volume fraction 
for the purpose of estimating the monomer fraction within our 
experimental data. The monomer fraction calculated using this 
power law for the experimental results is listed in Table 1 and is 
shown as the open circles in Figure 4. 

3.1. The C-H Stretching Region (3200-2900 cm- 1 ) 

Figure 1 shows the 3200-2900 cm -1 region of the spectrum 
of naphthalene in argon matrices at different concentrations and 
that of pure naphthalene. This region of the spectrum is dom- 
inated by two strong bands, both assignable to C-H stretching 


vibrations based on good agreement with the frequencies pre- 
dicted by DFT calculations using the 4-3 1G expansion basis set 
and B3LYP functional (Langhoff 1996). The spectral changes 
in going from the top to the bottom spectrum are caused by 
increasing interactions between naphthalene molecules as their 
concentration increases. These interactions perturb the vibra- 
tional force field, inducing changes in fundamental vibrational 
frequencies. 

The 3079 cm -1 band in the Ar/Nap = 1000 spectrum cor- 
responds to a vibration with b 2 U symmetry. With increasing 
fractional concentration of naphthalene, the b 2 t , C-H stretching 
mode broadens and is redshifted to 3068 cm -1 in pure naphtha- 
lene. The first five spectra in Figure 1 show that this 1 1 cm” 1 shift 
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Figure 4. Black squares are the simulated monomer fraction as a function of 
1 minus the volume fraction for initial intensities of 10 -6 , 5 x 10 -6 , 10 -5 , 
5 x 10 -5 , 10 -4 , 2.5 x 10 -4 , 5 x 10 -4 , and 10~ 3 A -3 . Statistical error bars are 
the size of the squares. The solid line is an unweighted linear fit to the data points 
constrained to pass through the origin and the fitted slope is 26.3 ± 0.5. The 
open circles are the monomer fractions calculated for the experimental results 
using this power law (see the text). 

(A color version of this figure is available in the online journal.) 


is not gradual, but rather caused by the growth of perhaps two 
distinct new features which produce a red shoulder that grows 
as the Ar/Nap concentration decreases. These bands, centered 
near 3075 cm -1 , overlap and ultimately overtake the narrower 
3079 cm” 1 matrix-isolated naphthalene band, which appears to 
be completely replaced when Ar/Nap = 6. It is significant that 
the peak position of the b 2 U C-H stretching band produced by 
naphthalene agglomerations at Ar/Nap = 6 (3074 cm" 1 ) is dif- 
ferent from that of pure naphthalene (3068 cm" 1 ) as shown in the 
bottom spectrum. This indicates that the structure of the naph- 
thalene agglomeration produced at very high concentrations and 
the corresponding intermolecular interactions are different from 
those in pure naphthalene. 

The 3066 cm" 1 band in the Ar/Nap = 1000 spectrum 
corresponds to a vibration with b\ u symmetry. The b\ u C-H 
stretching mode shows similar behavior to its nearby partner. 
The shoulder that grows with increasing concentration of 
naphthalene peaks at about 3062 cm" 1 , some 4 cm" 1 to the 
red of the original, and the same as the 4 cm” 1 offset of the 
red shoulder growing adjacent to the 3079 cm" 1 b 2 t , band. 
However, the 15 cm" 1 redshift between the b\ u vibration for 
the highly dilute spectrum versus that for pure naphthalene is 
some 5 cm" 1 greater than the 1 0 cm” 1 redshift for the 3079 cm” 1 
b 2 u vibration. 

The effect of increasing naphthalene concentration in the 
argon matrix can produce spectral shifts as follows. In the 
most dilute case (Ar/Nap = 1000, adjusted to 250 for this 
calculation), the fraction of naphthalene molecules that are 
isolated in the matrix is 0.90 ± 0.03. Consequently, since 
the infrared absorption bands arise primarily from isolated 
naphthalene molecules, the bands in the Ar/Nap = 1000 
spectrum are well defined and have widths on the order of 
a few cm” 1 . At Ar/Nap = 530, the monomer fraction drops 
slightly (0.82 ± 0.03), consistent with little or no change in 
the spectrum. The Ar/Nap = 200 spectrum shows a strong 
red wing on both features. In this case, the monomer fraction 
has dropped to 0.59 ± 0.08. At Ar/Nap = 160, this fraction 
becomes 0.51 ± 0.09, consistent with larger red shoulders. 


At Ar/Nap = 47, the monomer fraction is only 0.10 ± 0.04. 
Not surprisingly, it is at this concentration that the new bands 
producing the shoulders overtake the initial narrow features. 
At Ar/Nap = 6, our equation for the monomer fraction breaks 
down, since the volume fraction of naphthalene exceeds the 
maximal limit of 0.125 for the Mansson & Rudemo (2002) 
global thinning process. Ar/Nap = 6 clearly corresponds to 
massive agglomeration of naphthalene molecules within the 
matrix with a significant change in the spectrum in going from 
Ar/Nap = 47 to Ar/Nap = 6. There is no longer a contribution at 
the initial frequencies, but now both bands are symmetric with 
Gaussian profiles. There still remains a significant difference 
between the spectra of these agglomerations and that of pure 
naphthalene as shown by comparing the bottom two spectra in 
Figure 1. 

3.2. The C-H Out-of-plane Bending Region (900-400 cm -1 ) 

Figure 2 shows the 900-400 cm" 1 region of the argon matrix- 
isolated spectra discussed above along with the correspond- 
ing spectrum of pure naphthalene. For aromatic molecules, 
the strong bands in this region originate in C-H out-of-plane 
(CH oop ) bending vibrations and they are characteristic of the 
number of adjacent hydrogen atoms on each ring (e.g., Silver- 
stein & Bassler 1967; Bellamy 1960). A feature near 790 cm” 1 
dominates this region for naphthalene because it has only 
four adjacent (quartet) hydrogens per ring (e.g., Hudgins & 
Allamandola 1999). This is readily assigned to the b\ u CH oop 
bend, again based on good agreement with the frequency 
788.4 cm" 1 predicted by DFT calculations (Langhoff 1996). 
As with the C-H stretching bands discussed above, the spectral 
changes in going from the top to the bottom spectrum are caused 
by increasing interactions between naphthalene molecules as 
their concentration increases. However, for the CH oop bending 
mode, the band shift upon clustering is to the blue not the red 
as is the case for the C-H stretching modes. 

For the most dilute cases (Ar/Nap = 1000, 530, and 200), 
a strong, narrow band (784 cm" 1 , FWHM 0.9 cm" 1 ) accom- 
panied by one or more well-defined satellite bands dominates 
the entire region covered in Figure 2. The 784 cm” 1 band is 
readily assigned to isolated naphthalene and falls within 2 cm" 1 
of the position measured for gas-phase naphthalene near room 
temperature (Pirali et al. 2009). For the CH oop bending feature, 
the overlapping satellite bands overtake the isolated naphtha- 
lene 784 cm" 1 band at Ar/Nap = 160. This is not the case 
for the C-H stretching features until Ar/Nap = 47 (Figure 1). 
As with the band profile, width and position in the CH stretch- 
ing region, the spectrum changes significantly in going from 
Ar/Nap = 47 to Ar/Nap = 6. The narrow band and broader 
blue satellite feature are now replaced by a strong, Gaussian- 
like band peaking at 786 cm” 1 arising from agglomerations of 
naphthalene molecules. Interestingly, the peak position of this 
broad band is not significantly shifted from the central position 
of the narrow band and satellite features associated with the 
more dilute samples nor is there much difference with the band 
produced by pure naphthalene which peaks at 785 cm" 1 . 

Two other bands fall in this spectral window, one near 
620 cm" 1 and the other near 479 cm” 1 . Again, for the most 
dilute mixtures, these bands fall within 2 cm" 1 of those for 
gas-phase naphthalene (Pirali et al. 2009). Consistent with their 
origin in C-C-C skeletal vibrations, both are quite weak. Apart 
from slight broadening and at most a few cm" 1 shift, both 
of these bands show little change with increasing naphthalene 
concentration. 
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3.3. The C-C Stretching and C-H In-plane Bending 
Region ( 1550-900 cnT" 1 ) 

Figure 3 shows the 1550-900 cm -1 region of these matrix- 
isolated naphthalene spectra, a region rich in bands. For aromatic 
molecules, features between roughly 1700 and 1500 cm -1 arise 
from C-C stretching vibrations, those in the 1500-1250 cm -1 
region are associated with vibrations in which C-C stretching 
and C-H in-plane (CH; p ) bending motions are coupled, while 
those between 1250 and 1100 cm -1 primarily arise from un- 
coupled CHjp bending vibrations. The nature of these vibrations 
is again readily assigned by comparison with the frequencies 
predicted by DFT calculations (Langhoff 1996). The peak fre- 
quencies for the principle bands we measure for matrix-isolated 
naphthalene all lie within 5 cm -1 of their corresponding gas- 
phase counterparts as reported by Pirali et al. (2009). As above, 
the changes in going from the top to the bottom spectrum are 
due to increasing interactions between naphthalene molecules 
as their concentration increases. 

The discussions of the overall changes in the C-H stretching 
and bending regions with increasing naphthalene concentration 
discussed above can be applied to most of the changes shown 
in Figure 3. The bands are sharp and narrow in the most dilute 
case and broaden as satellite features grow in as the naphthalene 
concentration increases. With one exception, all of the naphtha- 
lene cluster bands fall, at most, a few wavenumbers to the red of 
the band associated with matrix-isolated naphthalene. While the 
bands broaden upon clustering. Table 2 shows that, for the ma- 
jority of these features, the bandwidths are consistently smaller 
than those for the C-H stretching and CH oop bands. The entries 
in Table 2 also show that the band shifts for most of the features 
in this region are modest, on the order of 1 or 2 cm -1 , similar to 
the shifts that occur with the CH oop modes. Taken together with 
the knowledge that most of the bands in this region arise from 
C-C stretching vibrations coupled with CH lp bending vibra- 
tions, the narrow bands and small shifts for most of the features 
in this region are a clear indication that clustering does not sig- 
nificantly influence the bonding within the naphthalene carbon 
skeleton nor does it influence the amplitudes of the different 
carbon-carbon vibrations that produce these bands. 

There is however a portion of this region of the spectrum 
which does show unusual changes as clustering increases. A 
narrow band (FWHM < 1 cm” 1 ) at 1131 cm” 1 is the only 
feature present between 1140 and 1100 cm” 1 in the spectrum 
of the Ar/Nap = 1000 sample. A blue wing grows in the 
Ar/Nap = 530 and 200 matrices. This starts to change con- 
siderably at Ar/Nap = 160. By Ar/Nap = 47, a red shoulder is 
clearly present, giving the impression that the sharp 1131 cm” 1 
and 1139 cm” 1 bands are now atop a broad underlying compo- 
nent. This behavior continues through the Ar/Nap = 6 sample, 
where the broad component spans the 1150-1120 cm” 1 range 
in a manner similar to that in the spectrum of pure naphthalene. 
Recall that these bands fall squarely in the CHi p bending region. 
Thus, the greatest changes in the 1550-900 cm” 1 region can 
be understood as arising from the stronger effect naphthalene 
clustering has on the peripheral C-H bonds than on the skeletal 
C-C bonds. 

4. ASTROPHYSICAL IMPLICATIONS 

It is important to keep in mind that this study is limited to 
clusters of neutral naphthalene, a very small aromatic molecule 
by interstellar and solar system standards, and care is warranted 
in applying these results to observational data. However, in view 


of the paucity of information on the spectroscopic effects that 
clustering has on PAH infrared spectra, a few interesting points 
emerge which merit some discussion in an astrophysical context. 

The CH stretching region. Since emission in the CH stretching 
region is dominated by the smallest members of the astronomical 
PAH population (Schutte et al. 1993), small PAH clusters 
containing on the order of 30-50 carbon atoms may contribute 
to the emission in this region. 

The CHi p bending region. A broad plateau that underlies 
the strong 7.7 /xm emission band complex (Cohen et al. 1986, 
1989) shows different spatial behavior than that of the narrower 
PAH emission features (Bregman et al. 1989). Rapacioli et al. 
(2005) carried out a detailed study of the spatial behavior of the 
emission in this region in two extended objects, NGC 7023 and 
p Oph-SR3. Focusing on the region near 7.7 /xm, they analyze 
the emission in terms of three components, neutral PAHs, PAH 
cations, and Very Small Grains (VSGs). They propose that PAH 
molecules are trapped in VSGs, perhaps as PAH clusters. 

In view of this, it is interesting that the strongest effect 
clustering has on the 1500-900 cm” 1 (667-11 /xm) region in 
the spectrum of naphthalene is an overall broadening of the 
features with red and blue wings stretching from 1150 to 
1100 cm” 1 (8.7-9 /xm), the region of the CH ip vibrations for 
quartet hydrogens. Quartet hydrogens are far less important 
than solo hydrogens in the astronomical PAH population. If this 
observed strong broadening for naphthalene clusters generalizes 
to PAH clusters with primarily solo hydrogens, one could 
speculate that clustering of interstellar PAHs might in general 
produce a series of broad features near the 7.7 /xm interstellar 
band since that is where the CHi p bending vibrations of solo 
hydrogens are most strong. 

The CH oop bending region. Two of the puzzles remaining 
in deciphering the astronomical spectrum in this region are 
understanding the peculiar blue wing of the astronomical 
787 cm” 1 (12.7 /xm) band (Hony et al. 2001; Bauschlicher 
et al. 2008) and the very broad plateau which underlies the 
entire region between roughly 900 and 700 cm” 1 (1 1-14.3 /xm) 
(Roche et al. 1989). 

The 12.7 /xm astronomical band is attributed to the CH oop 
bending vibration of astronomical PAHs with duo and trio 
hydrogens (see discussion in Bauschlicher et al. (2008), and 
references therein). However, whereas redshifting and band 
broadening are not unusual, blue shading is enigmatic among 
the astronomical PAH emission features. The observed blue 
shading of the naphthalene CH oop band to 786 cm” 1 (12.7 /xm) 
with decreasing Ar/Nap ratio is temptingly consistent with the 
phenomenon. 

If PAH clusters in general contribute to the blue component, 
one would predict spatial variation between the main peak and 
the blue wing. Since there seem to be few studies of the spatial 
behavior of the 12.7 /xm feature in extended sources, the work 
of Roche et al. (1989) on NGC 7027 and Rapacioli et al. 
(2005) on NGC 7023 is particularly important. Both studies 
indicate that the 12.7 /xm band and its blue wing show different 
spatial behavior with position in the object, implying that these 
two components track different members of the interstellar 
PAH family. Turning to the underlying plateau, Roche et al. 
(1989) clearly show that it tracks a separate component of the 
astronomical PAH family, one which extends well beyond the 
region in which the PAH features are strong. This behavior is 
similar to that reported by Rapacioli et al. (2005) for the bands 
in the NGC 7023. Clearly, the spectroscopic effects observed 
here must be tested with further studies on clusters involving 
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larger PAHs with different edge structures, mixed PAH clusters, 
and PAH clusters with different electrical charge, all species 
which likely make up the astronomical PAH population. Since 
PAH clusters are model components of Jupiter (Friedson et al. 
2002) and Titan’s atmospheres (Wilson & Atreya 2003), the 
information presented here may also be applicable to the 
spectroscopy of these objects. 

5. CONCLUSIONS 

This paper describes the changes in the infrared spectrum 
of naphthalene that occur as isolated naphthalene molecules 
agglomerate into clusters. As naphthalene clusters grow, bands 
shift and broaden. The strongest perturbations are associated 
with C-H vibrations, not skeletal C-C vibrations. Spectral 
shifts are modest, some 4-12 cm -1 , and bands are broadened 
to 5-10 cm” 1 . The possible contribution of homogeneous 
naphthalene clusters to the widespread interstellar infrared 
emission spectrum is discussed. Specifically, the likelihood that 
PAH clusters contribute to the emission plateau that underlies 
the 7.7 /Tin and 8.6 ji m features and the blue shading of the 
12.7 /xm emission band is considered. 
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